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the  sonic  line  and  the  given  profile  contour  an  exact  supersonic 
flow  field,  and  in  addition,  still  from  the  supersonic  region,  the 
changes  which  occur  in  this  field  at  the  current  sonic  line  and 
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numerically  the  kernel  of  an  integral  which  expresses  the  changes 
of  the  flow  field  in  terms  of  the  changes  of  the  shape  of  the  sonic 
line.)  This  information  is  used  to  formulate  a  boundary  condition 
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of  the  shape  of  these  two  lines.  With  these  boundary  conditions 
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elimination.  This  procedure  may  be  visualized  as  a  Newton-Raphson 
iteration.  The  report  discusses  this  procedure  in  general  terms 
and  derives  specific  formulae  for  the  conditions  at  the  sonic  line 
and  at  the  shock.  The  kernel  mentioned  above  will  have  singulari¬ 
ties.  The  nature  of  these  singularities  is  explored. 
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SECTION  I 
INTRODUCTION 


In  the  Murman  Cole  procedure  the  supersonic  and  the 
subsonic  regions  are  described  in  the  same  manner.  However, 
the  difference  stars  for  the  subsonic  and  the  supersonic 
regions  are  not  the  same  and  one  needs  special  measures  for 
the  treatment  of  the  sonic  points  and  of  the  shock  points. 

The  strength  of  the  Murman-Cole  approach  lies  in  the  fact 
that  these  different  requirements  are  combined  into  a  relatively 
simple  procedure  with  a  minimum  of  logical  decisions.  Some 
of  the  simplicity  is  lost  if  one  formulates  the  problems  in 
terms  of  finite  elements.  It  is  still  possible  to  retain 
the  same  elements  in  the  subsonic  and  in  the  supersonic  regions. 
The  differences  in  the  numerical  approach  for  the  subsonic 
and  for  the  supersonic  regions  are  brought  about  by  the  choice 
of  weight  functions.  This  is  probably  somewhat  more  compli¬ 
cated  than  the  difference  star  modifications  of  the  Murman- 
Cole  approach.  Further  complications  arise  when  one  expresses 
the  conditions  at  the  sonic  line  and  at  the  shock.  The  finite 
element  procedure,  in  order  to  be  economical,  must  use  elements 
which  are  larger  than  one  mesh  of  the  finite  difference  grid 
system.  In  the  Murman-Cole  procedure  the  location  of  the 
sonic  line  and  of  the  shock  is  only  roughly  defined.  To 
leave  these  curves  uncertain  within  the  finite  elements  would 
impair  the  accuracy.  A  procedure  of  this  kind  has  been  studied 
elsewhere  (Ref.  1) . 

One  may  ask  whether  it  is  really  advantageous  to 
use  the  same  description  for  the  subsonic  and  in  the  supersonic 
regions.  The  present  report  gives  the  basis  for  a  procedure 
in  which  the  supersonic  and  the  subsonic  regions  are  treated 
separately  with  the  provision  that  they  match  at  the  sonic 
line  and  at  the  shock.  At  the  sonic  line  one  requires  that 
the  potential  and  its  conormal  derivative  agree.  At  the 
shock  one  imposes  the  shock  conditions,  namely  continuity  of 
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potential  and  mass  flow  through  the  shock.  In  the  present 
study  we  shall  show  how  this  can  be  accomplished  by  a  Newton- 
Raphson  procedure,  in  which  one  computes  corrections  to  an 
available  approximation  (presumably  the  result  of  a  previous 
iteration  step) .  Convergence  can  be  expected  if  the  starting 
approximation  is  sufficiently  close.  In  the  supersonic  region 
the  flow  can  be  computed  rather  rapidly  by  a  marching 
procedure.  It  is  then  possible  to  evaluate  not  only  the 
basic  supersonic  flow  field  but  also  variations  to  it,  and  to 
derive  in  this  manner  boundary  conditions  for  the  subsonic 
region  along  the  common  boundary  of  the  subsonic  and  the 
supersonic  regions,  namely  the  sonic  line  and  the  shock  of 
the  (current)  approximation.  The  resulting  equations  for  the 
subsonic  region  (which  are  linear  because  one  carries  out  a 
Newton  Raphson  procedure)  are  then  solved  by  direct  elimination. 
Iterations  are,  of  course,  necessary  because  of  the  nonlinearity 
of  the  problem.  The  iterations  of  the  Murman-Cole  procedure 
serve  a  dual  purpose:  they  determine  corrections  to  the  flow 
field,  and  at  the  same  time,  because  of  immediate  updating, 
they  also  take  into  account  the  nonlinearities  of  the  problem. 

The  procedure  just  described  ought  to  converge  if  one 
has  an  approximation  which  is  sufficiently  close,  but  it 
requires  in  each  iteration  step  that  one  compute  the  super¬ 
sonic  region  a  number  of  times  for  different  assumed  shapes 
of  the  sonic  line.  One  might  think  of  a  different  approach 
in  which  the  supersonic  region  is  computed  only  once  in  each 
iteration  step.  One  starts  with  an  approximation  to  the 
subsonic  field  and  determines  the  location  of  the  sonic  line. 
With  this  result  one  would  compute  the  supersonic  field;  in 
particular,  the  values  of  the  conormal  derivative  of  the  per¬ 
turbation  potential  at  the  sonic  line.  For  these  values  and 
with  the  shock  conditions,  corrections  to  the  subsonic  region 
could  be  computed.  This  would  complete  one  iteration  step. 
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Here  one  does  not  deal  with  a  Newton  procedure.  The  con¬ 
vergence  of  such  iterations  is  uncertain  even  if  one  has  a 
close  approximation,  but  if  the  method  converges,  then  the 
matching  conditions  at  the  sonic  line  will  be  satisfied. 
Whether  one  has  convergence  and  whether  the  procedure  is 
economical  must  be  determined  by  numerical  experiment. 
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SECTION  II 
THE  SUBSONIC  REGION 


It  is  assumed  that  an  approximation  for  the  subsonic 
part  of  the  flow  field,  including  an  approximation  for  the 
location  of  the  sonic  line  and  of  the  shock,  is  available. 

The  sonic  line  and  that  part  of  the  shock  which  leads  to 
subsonic  velocities  constitute  the  common  boundary  between 
the  subsonic  region  and  the  supersonic  region.  We  refer  to 
these  data  as  current  approximations.  It  is  assumed  that 
corrections  to  the  subsonic  part  of  the  flow  field  are 
obtained  by  a  finite  element  procedure.  To  fix  the  ideas 
we  consider  rectangular  elements.  For  the  present  discussions 
which  are  concerned  with  the  basic  questions,  the  character 
of  the  shape  functions  to  be  chosen  (bilinear,  biquadratic, 
or  bicubic)  is  unimportant.  The  current  sonic  line  and  the 
current  shock  cut  through  these  elements.  This  must,  of 
course,  be  taken  into  account  in  the  computations.  The 
subsonic  problem  can  be  formulated  by  means  of  Bateman's 
extremum  principle;  the  functional  whose  variation  must 
vanish  is  the  integral  of  the  pressure  over  the  region  under 
consideration.  A  linearization  is  necessary  because  the 
relation  between  the  pressure  and  the  velocity  gradient  is 
of  a  higher  degree  than  the  second.  The  extremum  formulation 
is  used  only  to  derive  the  governing  equations.  The  conditions 
at  the  shock  and  at  the  sonic  line  which  occur  here  are  not 
of  a  nature  which  allows  one  to  show  that  the  functional 
described  above  has  an  extremum.  Convergence  can  be  expected 
because  one  carries  out  a  Newton-Raphson  procedure  if  one 
has  a  sufficiently  close  starting  approximation. 

If  one  assumes,  temporarily,  that  along  the  current 
sonic  line  and  the  current  shock  the  conormal  derivative  to 
the  perturbation  potential  is  prescribed,  then  the  potential 
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along  this  line  appears  as  unknown;  the  conormal  derivative 
appears  in  the  equations  obtained  from  the  application  of  the 
extremum  principle  (which  would  then  hold)  in  the  inhomogeneous 
part.  (The  other  contributions  to  the  inhomogeneous  part 
arise  because  in  general  the  current  approximation  does  not 
satisfy  the  original  nonlinear  equations  exactly.)  Actually,  the 
conormal  derivative  is  unknown.  Computations  to  be  carried 
out  in  the  supersonic  region  express  it  and  also  the  change 
of  the  potential  along  the  sonic  line  and  along  the  shock  in 
terms  of  the  change  of  the  shape  of  the  sonic  line  and  of  the 
shock.  With  these  additional  relations,  one  computes 
simultaneously  the  change  of  the  potential  in  the  subsonic 
part  of  the  flow  field,  the  change  of  the  shape  of  the  sonic 
line,  and  the  change  of  the  shape  of  the  shock.  The  character 
of  the  conditions  at  the  sonic  line  and  at  the  shock  and  the 
computations  necessary  in  the  supersonic  region  will  be 
discussed  in  the  next  sections. 
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SECTION  III 

CONDITIONS  AT  THE  SONIC  LINE 

It  is  our  goal  to  express  the  change  of  the  potential 
and  of  its  gradient  in  terms  of  some  characterization  of 
the  deformation  of  the  sonic  line.  We  confine  our  attention 
to  the  supersonic  region.  Consider,  simultaneously,  a  basic 
flow  field  and  this  basic  flow  field  with  a  superimposed 
perturbation.  Points  for  which  the  velocity  vector  is  the 
same  in  the  two  fields  are  called  corresponding  points.  The 
vector  leading  from  a  certain  point  of  the  basic 
flow  field  to  its  corresponding  point  in  the  perturbed  field 
is  called  the  displacement  vector  of  this  point.  It  is 
practical  to  express  perturbations  at  the  current  sonic  line 
in  terms  of  the  component  of  the  displacement  vector  normal 
to  the  streamlines.  (For  the  simplified  transonic  equation 
this  is  the  component  in  the  y  direction.)  Let  s  be  a 
parameter  which  changes  monotonically  along  the  sonic  line 
(in  most  cases  the  x  coordinate  can  serve  as  such  a  parameter) 
and  let,  for  some  neighboring  field,  dn(s)  and  dfc(s)  be  the 
components  of  the  displacement  vector  for  points  of  the  sonic 
line,  respectively,  in  the  direction  normal  and  tangential 
to  the  streamlines.  The  computations  in  the  supersonic 
region  serve  to  express  dt(s)  in  terms  of  dn(s).  This  will 
be  discussed  later.  We  express  the  change  of  the  gradient 
of  the  potential  in  terms  of  dn(s)  and  dfc(s)  .  First  we 
compute  some  derivatives  pertaining  to  the  basic  flow  field 
in  terms  of  geometric  quantities.  This  can  be  done  without 
difficulty  even  for  the  full  potential  equation  in  a  local 
coordinate  system  where  the  temporary  x  axis  coincides  with 
the  flow  direction  and  the  y  axis  is  normal  to  it.  Then  one 
has  <f>  =  0  and  <j>x  =  a*  (because  we  consider  a  point  of  the 
sonic  line) .  The  potential  equation  yields 

4>  =  0  (1) 

yy 
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Let  R  be  the  radius  of  curvature  of  the  streamline  at  the 
point  in  question.  Then  one  has 

<+>  =  a*  •  -  (2) 
vxy  R 

a*  is  the  sonic  velocity,  usually  one  chooses  a*  =  1.  Carrying 
a*  within  the  equations  facilitates  a  check  of  the  dimensions. 
Let  6  be  the  angle  of  the  sonic  line  with  the  local  x  axis. 

At  the  sonic  line  one  has 

2  2  2 
4>x  +  V  =  a* 

and  since  in  the  chosen  system  of  coordinates  cj)^  =  0  one 
obtains  by  differentiation  with  respect  to  x  (along  the  sonic 
line) 

*x(<t>xx  +<»xy  tg  B)  =  0 


Hence 


*xx  =  ^xy  t(*  6  =  -a*R 


tg  B 


In  the  local  system  of  coordinates  dt(s)  and  d^(s)  are  the 
displacement  components,  respectively,  in  the  x  and  y 
directions.  The  components  of  the  gradient  of  the  potential 
in  the  directions  tangential  and  normal  to  the  stream  lines 
are  given  by 


Let  the  perturbation  potential  be  denoted  by  A<f>.  Expressing 
the  fact  that  at  corresponding  points  of  the  sonic  line  the 
velocity  vectors  in  the  basic  and  in  the  changed  flow  field 
are  the  same,  and  introducing  a  linearization  one  obtains 


A<j).  +  <J>  d  +  <j>  d  =  0 
Yt  Txx  t  Yxy  n 


Ad>  +  d>  d.  +  <f>  d  =  0 
Yn  Yxy  t  Yyy  n 


7 


Hence 


A4»t (s)  =  a*R-1  (tg  3  dfc(s)  -  dn(s))  (4) 

A*n(s)  =  -a*R_1dt(s)  (5) 

These  equations  constitute  boundary  conditions  for  the 
perturbations  in  the  subsonic  region  at  the  current  sonic 
line  provided  that  dt<s)  is  known  as  a  function  of  dn(s). 
This  relation  will  be  written  in  the  form 

s 

d.  (s)  =  /  K (s , t) d ' ( t) dt  (6) 

t  Q  n 

The  derivative  d  '  =  d(d  (t))/dt  is  used  in  order  to  obtain 
n  n 

a  kernel  K(s,t)  which  is  less  singular  (and  therefore  better 
suited  for  a  numerical  evaluation) .  The  determination  of 
K ( s , t)  will  be  shown  in  Section  V. 


SECTION  IV 

CONDITIONS  AT  A  SHOCK 


It  will  be  assumed  in  the  coming  dicussions  that  the 

shock  is  weak.  Then  the  change  of  the  entropy  through  the 

shock  is  small,  and  the  flow  downstream  of  the  shock  is  again 

a  potential  flow.  At  the  shock  the  potential  and  the  mass 

flow  through  the  shock  are  continuous.  (Continuity  of  the 

potential  assures  that  the  velocity  component  tangential 

to  the  shock  is  continuous.)  Quantities  referring  to  the 

flow  field  upstream  or  downstream  of  the  shock  will  have, 

respectively,  superscript  1  and  2.  It  is  assumed  that  at 

the  current  shock  position  the  basic  flow  field  is  known 

and  that  A<{>^,  A<}>  ^  and  A<J>  ^  are  expressed  as  functions 

x  y 

of  d  (s) ;  so  far  d  (s)  is  unknown.  The  basic  flow  field  in 
n  n 

the  supersonic  region  will  be  determined  for  the  current 
shape  of  the  sonic  line  by  means  of  the  method  of  character¬ 
istics.  The  potential  at  the  current  shock  position  upstream 

and  downstream  of  the  shock  are  usually  not  the  same  because 
(2) 

the  potential  ^  is  given  by  the  current  approximation  to 
the  subsonic  part  of  the  flow  field,  while  the  potential 
<J>  ^  is  taken  from  the  newly  computed  supersonic  region. 

The  shock  position  in  the  original  and  in  the  perturbed  flow 
fields  are  characterized  by 

x  *  ’‘shock1*1  and  *  =  xshock(y)  +  4xshockly) 


Expressing  the  requirement  that  at  the  changed  shock  position 
the  potential  is  continuous,  one  obtains 


'»(1)(xshock'yl 


^U)  lxshock'y) 


(1) 


^x  ^xshock ' Axshock ^ 


*(2)<x 


shock ' 


y)  +  6*121 (xshock'y) 


(2) 


*x  (xshock'y)Axshock(y) 
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Hence 


..  (2)  =  (1)  . 

^ shock  ^  xshock 


,y)  -  <t> (2)  (X 


shock 


,y)  +  l*shock.y) 


(**l) (xshcck'yl  ‘  (xsh°ck'y))Ax  ly)  shock  171 


The  only  unknown  quantities  in  the  right  hand  side  are  A 4>  ^ 
which  depends  upon  dn(s)  and  Axsj10ck^^'  For  simPlicity  t^ie 
equations  of  conservation  of  mass  for  the  flow  through  the 
shock  are  derived  only  for  the  simplified  transonic  equation; 
it  expresses  conservation  of  mass  in  terms  of  a  perturbation 
potential  which  characterizes  the  deviation  from  a  parallel 
flow  with  the  sonic  velocity.  This  differential  equation 
is  given  by 


-  (y  +  1)  a 


*-l 


X  XX 


yy 


=  o 


(8) 


Hence 


//Hy  +  U  a*-1<j>x4>xx  +  4>yy]dxdy  =  0 
and  by  partially  carrying  out  the  integrations 

f  -  [  (y  +  1 )  / 2  a  _14>2  dy  +  (f>  dx]  =  0 

9ft  x  y 

Hence  for  the  deviations  of  the  mass  flow  from  a  parallel 
sonic  flow  for  a  line  element  in  the  flow  field  with  slope 
dx/dy 

-  [  (Y  +  l)/2  a  +  ip  dx/dy]  dy  (9) 

x  y 

The  requirement  that  this  deviation  be  the  same  upstream  and 
downstream  of  the  shock  gives 
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h  +  l)a*"V2)  A*<2)  +  A<t(2)(dx 


shock 


/dy) 


<  <  .  *-l-.  (1)2  ,(2)2  (1) 

(  (y  +  l)/2)a  (4>x  -  <*x  )  +  (<t>y 


(2) 


(d*shock/dy> 


+  (y  ♦  lla*-1*'11!*™  4-  4*'1’  (dxshock/dy) 


A  *-l,.(l)  .  (1)  . ( 2)  ( 2)  .  . 

+  ( y  +  i ) a  (4>x  $xx  -  4>x  4>xx  )Axshock 


+  (<t> 


(1)  - 
xy 


Axshock(dxshock/dy) 


+  -  tfy2)  )  (dAXShOCk/dy) 


(10) 


The  unknown  quantities  on  the  left  hand  side  are  A <p  (2)  and 

(2)  x 
Atf>  .  (Actually  the  expression  on  the  left  is  the  conormal 

Y  ( 2 ) 

derivative  pertaining  to  the  differential  equation  for  A$  )  • 

The  unknown  quantities  on  the  right  hand  side  are  Axg(y)  and 

dn(s).  The  latter  function  occurs  in  A<J>  ^  and  . 
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ELk. 


SECTION  V 

TREATMENT  OF  THE  SUPERSONIC  REGION 


For  the  computations  of  the  supersonic  flow  field  we 
consider  the  shape  of  the  sonic  line  as  given.  In  addition 
one  knows  the  profile  contour.  In  this  section  we  shall 
discuss  the  application  of  the  method  of  characteristics, 
although  there  exist  finite  difference  methods  which  could 
also  be  used.  In  multidimensional  problems  finite  difference 
methods  would  probably  be  preferable. 

The  method  of  characteristics  for  the  full  potential 
equation  and  also  for  the  simplified  transonic  equation  is 
well  known.  To  express  the  ideas  which  are  important  in  the 
present  context,  we  shall  repeat  the  main  equations  for  the 
simplified  transonic  equation  but  without  reference  to  the 
underlying  mathematical  theory.  One  starts  from  Eq.  (8) , 
namely 

-(y  +  l)a  4>  +  d>  =0 

' 1  yx  xx  ^yy 

This  equation  is  rewritten  in  the  form 


[+(Y+l)1/2a*"1/24>^/2 


j^}  f+(Y+l) 1/2a*~1/2<2/3) 


For  curves  within  the  flow  field  with  the  slope 


dx 

dy 


+  (Y  +  D1/2a*-1/2^/2 


(11) 


this  equation  becomes 

[+(Y+l)1/2a*-1/2(2/3)<^/2  +  ♦y]  =  0  (12a) 


Hence,  for  these  curves 


const 


(12b) 


Curves  whose  slope  is  given  by  Eq.  (11)  are  called  characteristics . 

Eqs.  (12a)  and  (12b)  are  called  compatibility  conditions  for 

the  characteristics.  They  express  that  along  such  curves 

$  and  <}>  cannot  be  chosen  independently.  One  speaks  of  left 
x  y 

or  right  going  characteristics  to  suggest  their  orientation 
with  respect  to  the  streamlines.  (The  positive  sign  in  Eq.  (11) 
refers  to  left  going,  the  negative  sign  to  right  going 
characteristics.)  In  the  method  of  characteristics  one  computes 
a  net  formed  of  such  curves.  The  state  at  a  grid  point  of 
such  a  net  is  readily  found.  Assume  that  one  knows  <J>x  and  <J> 
at  the  points  1  and  2  of  Figure  1.  Points  1  and  2  must  not 
lie  on  the  same  characteristics.  Let  point  3  be  the  point 
where  two  characteristics  of  different  families  through  the 
points  1  and  2  intersect.  Then  one  has 


-<Y+1)  1/2a*‘1/2(2/3H^  + 


...1/2  *-1/2,  ,  ...  .3/2  .  . 

+  (y  +  l)  a  (2/3)<J>x'3  +  <(>  3 


-  (y+1)  1/2a*-1/2  (2/3)  +  <t>  , 

A  t  L  y  i  x 

(13) 

.  ,  ...1/2  *-1/2.-.,,  .3/2 
+  (y+l)  a  (2/3)4»Xf2  +  $Vi2 


(Second  subscripts  1,  2  or  3  refer  to  the  numbering  of  the 
points.)  The  right  hand  sides  are  known.  Therefore,  one 
can  compute  immediately,  <{>x  3  and  <J>  3-  After  these  quantities 

have  been  determined,  one  finds  x3  and  y3,  by  applying 
Eqs.  (11)  in  difference  form 


x 

y 

x 

y 


3 

3 

3 

3 


(Y+1)  1/2a*"1/2(l/2)  (4-J-7?  +  $x7?) 


-(Y+l)1/2a*'1/2(l/2)  <♦*^3  +  ^21 


(14) 


The  flow  field  in  the  vicinity  of  the  sonic  line  is 
regular;  that  is,  4>  can  be  developed  in  power  series  with 
respect  to  x  and  y  (unless  a  singularity  propagating  along  a 


13 


characteristic  occurs  at  the  point  for  which  the  development 

is  to  be  carried  out) .  However,  the  shape  of  the  characteris- 

1/2 

tics  is  singular,  because  of  the  power  <f>  '  in  Eq.  (11)  . 

This  can  be  taken  into  account  when  one  forms  Eqs.  (14)  in 
the  segment  of  the  characteristics  next  to  the  sonic  line. 

If  one  chooses  a  point  of  the  sonic  line  as  origin  of  an 
x,  y  system  and  aligns  the  x  axis  with  the  local  flow 
direction,  then  the  lowest  order  term  in  the  development  of 
<(>  is  given  by 


hence 


=  const  xy 


<t>  =  const  y 


Therefore,  for  the  characteristics  in  this  vicinity 

dx  r,  . , , 1/2  *-1/2  1/2  1/2 

=  +  (  y  +  1  )  a  const  y 

and  for  the  characteristics  through  the  origin 

-.  .  .  1/2  *-1/2  1/2  .  _  3/2 

x  =  + ( y+1 )  a  '  const  (2/3)y 


Expressing  y  by  $  ,  one  obtains 


— ,  .  .  1/2  *-1/2  4.-1/1/H*  3/2 

x  =  +(y+l)  a  const  (2/3)<fx 


Hence,  for  two  points  of  this  characteristic 

*3/2  -  »3/2 

X1  X2  -  1/2  *-1/2  X1  X2 
"  *  +(y  +  l)i7  a  X/ ^  ( 2/3 )  —  7 


Yl  ’  y2 


- 


(15) 


1/2 

Eqs.  (14)  and  (15)  are  exact  if,  respectively,  and  $ 

and  linear  functions  of  y.  Eq.  (15)  therefore,  gives  an 
acceptable  approximation  for  the  secant  directions  of  the 
characteristics  throuqhout  the  flow  field,  while  Eqs.  (14) 
becomes  inaccurate  in  a  neighborhood  of  the  sonic  line. 
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The  sonic  line  is  reached  by  only  one  downstream  going 

characteristic  (Fig.  2) .  The  state  at  the  sonic  line  is 

determined  by  the  compatibility  condition  for  this  characteristic 

and  the  condition  that  at  the  sonic  line  <b  =0. 

x 

Only  one  downstream  going  characteristic  arrives  at 
the  surface  of  the  profile  (Fig.  3) .  For  such  a  point,  one 
has  the  condition  for  the  slope  of  this  characteristic,  the 


compatibility  condition,  and  as  the  boundary  condition,  the 


value  of 


at  the  contour . 


(y+1)  1/2a*-1/2  (2/3)  +  4>  2  =  (Y+D1/2a 

y  #  ^ 


*-1/2,-.-,  3/2  ^ 

(2/3)<2  +  *y,l 


.3/2  _  3/2 

-(1/2)  (y+1  )1/2a*~1/2  |  -Jiil - 

3  ^x , 2  "  *x,l 


*  2  =  ftXj)  (16) 

From  these  three  equations,  <4>  _ ,  <J>  -  and  x-  are  found  by  an 

x  t  £  y  f  £  ^ 

iteration.  To  carry  out  the  method  of  characteristics  one 
also  needs  data  along  an  initial  line.  They  are  given  by  an 
analytic  expression  valid  for  the  entrance  corner  of  the 
supersonic  region.  This  will  be  discussed  later.  Busemann 
makes  a  special  choice  of  the  values  of  and  <J>x  which  occur 
at  the  grid  points  of  the  characteristic  net.  This  allows 
him  to  work  with  predetermined  characteristic  slopes.  In 
the  present  context,  it  is  probably  preferable  to  forgo  this 
simplification.  One  must,  however,  make  sure  that  the 
differences  in  <p^  for  adjacent  points  of  the  characteristic 
grid  are  not  too  large  (because  of  accuracy  requirements) 
and  not  too  small  (in  order  to  avoid  unnecessary  work) .  The 
first  condition  will  be  the  reason  that  one  must  start  new 
downstream  going  characteristics  at  the  surface  of  the 
profile. 
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f 


To  obtain  the  local  changes  A<|>  and  A<f>  one  might  start 

x  y 

with  the  linearized  form  of  Eq.  (7) , 

-(Y+l)a*"14.xA(^xx  -  (Y+l)a*'1A<j>x<{.xx  +  A({.yy  =  0  (17) 

The  characteristic  conditions  are  obtained  by  rewriting  this 
equation  in  the  form 

+  [I(Y+l)1/2a*'1/2^/2  ^  -  <Y+l>a*-\x4<,x  -  0 

Here  all  quantities  except  At|)x  and  A 4> y  are  considered  as 
known.  The  characteristics  are  again  in  line  with  the  slope 


dx/dy  =  +(Y+l)1/2a*~1/2<i>  1/2 


that  is,  they  are  identical  with  the  characteristics  of  the 
original  field.  The  compatibility  conditions  assume  the 
form 

+  (Y+1)  1/2a*“1/2^1/2  A<(>  )  +  ^-(A <}>  )  -  (y+1)  a* —1  cf>  A(j)  =  0 

-  x  dy  Yx  dy  y  xx  x 


The  singularity  which  occurs  in  the  shape  of  the 

1/2 

characteristics  and  the  factor  <p  which  occurs  in  the 
last  equation  would  make  further  discussions  for  the  vicinity 
of  the  sonic  line  necessary. 

It  is  simpler  if  one  carries  out  the  method  of 
characteristics  for  the  neighboring  field  as  well  as  for  the 
original  field  without  introducing  a  linearization.  At  the 
beginning  of  Section  II  we  defined  the  displacement  of  a 
point  the  flow  field  due  to  the  superposition  of  a 
perturbation:  the  displacement  vector  is  the  difference  in 


the  local  vectors  of  points  of  the  perturbed  flow  field  and 

of  the  original  flow  field  for  which  the  state  is  the  same. 

In  applying  the  method  of  characteristics  for  the  neighboring 

field  it  is  best  to  compute  these  displacement  vectors.  The 

hodograph  equation  is  linear,  from  this  fact  it  follows 

that  one  obtains  again  an  exact  solution  in  the  physical 

plane,  even  for  the  full  potential  equation,  if  one  multiplies 

all  displacement  vectors  by  a  constant  factor.  As  far  as 

the  displacement  vectors  are  concerned,  this  makes  a  limiting 

process  to  an  infinitesimally  small  perturbation  unnecessary. 

A  limiting  process  is,  however,  necessary  to  determine  from 

the  displacement  vectors  the  local  changes  A(J>  and  A<J>  .  Let 

x  y 

Ax  and  Ay  be  the  components  of  the  displacement  vector. 

From  the  condition  that  at  points  of  the  original  and  of 
the  perturbed  field  which  correspond  to  each  other  the 
velocity  vector  is  the  same,  one  then  obtains 


The  derivatives  <{>  ,  <i  and  it  are  taken  from  the  basic 

xx  Yxy  yy 

flow  field.  The  perturbations  A<!>x  and  Atf^  are  needed  only 
at  the  current  location  of  the  sonic  line  and  at  the  current 
location  of  the  shock. 

The  method  of  characteristics  which  determines  primarily 
the  location  of  the  grid  points  for  a  given  state  of  the 
flow  gives  the  displacement  vectors  directly.  To  start  the 
procedure,  one  must  use  the  displacement  vectors  along  some 
initial  line  where  they  will  be  computed  analytically  from 
data  at  the  entrance  corner  of  the  supersonic  region. 

If  one  has  points  1  and  2  in  the  original  field  and 
points  1'  and  2'  in  the  neighboring  field,  then  one  obtains 
the  ooint  3'  in  the  neighboring  field  by  drawing  the  parallel 


to  the  characteristics  13  and  23  (see  Figure  4) .  The  displace¬ 
ment  vector  is  then  given  by  the  line  33'.  Actually  it  can 
be  computed  if  only  the  displacement  vectors  1,1'  and  2,2' 
and  the  slopes  of  the  characteristic  1,3  and  2,3  are  known. 

Let  point  2  be  the  point  of  the  original  flow  field 
where  the  downstream  going  characteristic  through  point  1 
of  the  perturbed  field  which  corresponds  to  point  1.  For 
the  displacement  vector  at  point  2  the  component  d^  normal 
to  the  streamlines  is  known.  Then  one  has  as  locus  for  point 
2'  the  parallel  through  point  1*  to  the  characteristic  1,2 
and  the  parallel  to  the  streamline  at  point  2  at  a  distance 
d^.  With  these  data  one  determines  the  component  of  the 
displacement  dfc  in  streamwise  direction  at  point  2.  For  this 
computation,  only  the  displacement  vector  1,1',  the  assumed 
displacement  component  d^  at  point  2  and  the  slope  of  the 
characteristic  1,2  are  needed. 

At  the  profile  contour  one  must  express  the  fact  that 
the  local  flow  angle  remains  unchanged  if  one  superimposes 
a  perturbation.  One  has,  in  general,  for  corresponding 
points 


A0  +  (d0/dn)dn  +  (d9/dt)dt  =  0 
Therefore,  at  the  profile  contour 

(d0/dn)dn  +  (d0/dt)dt  =  0  (20) 

d0/dt  is  the  local  curvature  of  the  profile.  d0/dn  is  computed 
from  the  basic  flow  field,  best  probably  from  the  data  along 
the  two  characteristics  which  pass  through  the  point  of  the 
contour  in  question.  Let,  temporarily,  the  y  direction 
coincide  with  the  normal  to  the  streamlines  (and  to  the 
profile  contour).  Then  one  has,  by  differentiation  along 
the  characteristics  1,2  and  1,3 


(de/dy)  =  de/dn  +  d©/dt (dt/dy) %  ^ 

where  d0/dy  and  dt/dy  is  to  be  formed  along  the  respective 
characteristics  (subscripts  £  or  r  refer  to  left  or  right 
going  characteristics) .  These  equations  are  formed  for  both 
characteristics  at  the  point  of  the  profile  contour.  The 
derivatives  (d0/dy) ^  and  (d0/dy)r  for  the  point  at  the  contour 
are  found  from  the  basic  flow  field.  Making  use  of  the  fact 
that  (dt/dy)  =  - (dt/dy)  ,  one  obtains 

d0/dn  =  1/2 ( (de/dy) ^  +  (d0/dy) r) )  (21) 

of  course,  d^  and  d  (in  Eq.  (20))  must  be  expressed  in  terms 
of  Ax  and  Ay.  One  has,  for  the  simplified  transonic  equation 

dn  =  Ay  ,  dt  =  Ax  (22) 


Point  2'  lies  on  the  parallel  to  the  characteristic  2  through 
point  1',  and  on  the  straight  line  defined  by  Eq.  (20). 

The  flow  field  in  the  entrance  corner  to  the  supersonic 

region  is  determined  analytically.  At  the  sonic  line,  one 

has  p  =0.  Let  the  sonic  line  be  given  by  y  =  y  (x)  . 

x  n  j  j  j sonic 

<p  +  <p  (dy  /dx)  =  0 

xx  Yxy  1  sonic 

Here,  <|>  is  given  by  the  boundary  condition  at  the  profile. 
In  geometric  terms  one  has 

*-l  -1 

a  <j>  =  R  (23) 

Yxy 


where  R  is  the  radius  of  curvature  at  the  point  under  con¬ 
sideration.  Further  information  is  provided  by  the  observa¬ 
tion  that  one  can  obtain  the  development  of  the  potential 
for  the  vicinity  of  a  point  of  the  contour  if  the  velocity 
vector  at  the  contour  is  given,  for  the  contour  is  not  a 
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characteristic  line.  In  the  present  context,  the  result 
is  immediately  obvious.  One  obtains  from  the  simplified 
transonic  equation  for  the  entrance  corner  of  the  supersonic 
region 


0 


One  then  obtains  in  a  sufficiently  small  neighborhood  a 
coordinate  system  whose  origin  lies  at  the  entrance  corner 
of  the  supersonic  region 


<(>  =  4>  (-(dy  .  /dx)x  +  y) 

Yx  Txy  ■'sonic'  1 


b  =  <f>  (0 , 0)  +  <p  x 
y  y  yxy 


(24) 


To  start  the  procedure  one  will  compute  the  value  of  <J> 
at  point  1  of  Figure  7;  4>x  =  0  because  this  is  a  point  of 
the  sonic  line. 

Next,  one  introduces  a  perturbation  by  assuming  that 
for  y  =  0  (that  is  along  the  contour)  4>x  is  given  by 


<)>  (x,0)  =  A<f>  (0,0)  +  4>  (0,0)x  +  A<t>  (0,0)x 

a  A  A  A  AA 


Since  ()>Xy  is  determined  by  the  profile  shape,  one  has  in  the 
changed  flow  field 


<l>x(x,y)  +  A<t>x(X'y)  =  A  4>x  (0,0)  +  <j>xx(0,0)x  +  A<f>xx(0,0)x  +  $  ( 

Thus,  for  the  local  change 

^4>x(x,y)  =  A4>x  (0,0)  +  A  c^xx  ( 0 , 0)  x  (25) 

Now,  we  invoke  the  differential  equation  (for  x  =  0,  y  =  0) 

-(Y+l)a*-1A<f >  (0,0)  [4  (0,0)  +  A*vv(0,0)]  +  A*  (0,0)  =  0 

x  xx  xx  yy 


» 0)  y 
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Hence,  with  a  linearization 


A<t>yy(0,0)  =  (y+l)a*"1A<{)x(0,0)(J)xx(0,0)  (26) 

and 

A<j>  (x,y)  =  (Y+l)a*_1A4)  (0,0)  <j>  ( 0 , 0)  y  (27) 

y  xjcx 

To  determine  from  these  results  the  displacement  vector 
(components  Ax  and  Ay)  in  the  flow  field,  we  express  the  fact 
that  at  corresponding  points  the  values  of  <J>x  and  <}>y  are  the 
same. 

<J>  (0,0)x  +  <t>  ( 0 , 0)  y  =  A6  (0,0)  +  4  (0,0)  (x+Ax)  +  A4  (0,0)x 

aa  x  y  x  xx  xx 


+  4>x  y(0,0)  (y  +  Ay) 


<t>xy(0,0)x 

Hence 
Ax (x,y)  = 


d>xy(0,0)  (x  +  Ax)  +  Ad>yy(0,0)y 

A4>  (0,0) 

t  (0T0Y  Y 

xy 


Ay (x,y) 


1 

< p  (0,0) 

vxy 


[-A0x (0,0) 


A4>xx(°f0)x  + 


4>xx  (0,0)  A  (0,0) 

^xy(0,0) 


y) 


(28) 


(29) 


Here  Ad)  (0,0)  is  expressed  by  Eq.  (26).  So  far  unknown  are 

yy 

the  assumed  perturbations  at  the  contour  A<|>  (0,0)  and  A4  (0,0). 

The  purpose  of  these  computations  is  the  evaluation  of  the 

kernel  K(s,t)  in  Eq.  (6).  One  chooses  dn(t)  =  0,  up  to  a 

point  t  =  t  and  d  (t)  =1,  from  there  on;  in  other  words, 

n 
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dn‘ (t)  and  <$(t  -  6).  For  this  choice  one  computes  dt(s,t). 
This  function  depends,  of  course,  upon  the  value  of  t .  Then, 
from  Eq.  (6) 


dt(s,x)  =  K ( s , x ) 


(30) 


These  computations  must  be  carried  out  for  a  sufficient  number 
of  values  of  t  and  s  so  that,  if  necessary,  intermediate 
values  of  K(s,t)  can  be  obtained  by  interpolation. 

If  the  point  t  =  t  lies  at  the  entrance  corner  of  the 
supersonic  region,  (x  =  0,  y  =  0  in  Eq.  (29))  then  one  obtains 

Ay ( 0 , 0)  =  -A$x<0,0)/>  (0,0) 


Hence  A*x(0,0)  =  “4>xy(0,0)  Ay(0,0) 


(31) 


For  the  evalution  of  K(s,t)  described  above,  one  has  d  (x)  = 

n 

Ay  =  1  also  for  other  points  of  the  sonic  line.  Therefore, 
from  Eq.  (29) 


1  = 


xy 


ra775T['aV0-01  -  <a*x*<°-0)  - 


<t>  (0,0)  Ad>  (0,0)  dy 

'xx  '  vyy  *  sonic 


<t>  (0,0) 

rxy 


dx 


)  x 


Hence,  with  Eqs.  (31)  and  (26) 

A4>xx(0'°>  *  “  ( Y+l)  a*-1^xx  ( 0 , 0)  2 


dy 


sonic 

dx 


Ay(0,0) 


This  allows  one  to  compute  Ax(x,y)  and  Ay(x,y) ,  for 
the  point  at  which  the  method  of  characteristics  starts 
(points  1  in  Figure  7)  by  means  of  Eqs.  (25)  and  (27). 

The  remaining  discussions  address  themselves  to  the 

character  of  the  kernel  K(s,t)  for  points  at  the  sonic  line, 

other  than  the  initial  point.  The  computation  described 

above  introduces  at  the  station  t  =  t  a  jump  of  d  (t) ,  one 

n 

therefore  expects  that  dfc(s)  will  be  singular  for  s  =  t. 
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It  is  desirable  to  know,  at  least  qualitatively,  the 
character  of  these  singularities,  for  the  numerical  evalua¬ 
tion  gives  only  pointwise  data.  For  this  purpose  it 
is  best  to  use  the  simplified  hodograph  equation,  which 
arises  from  the  simplified  transonic  equation  by  the  Legendre 
transformation.  To  make  the  report  self-contained,  the 
fundamental  derivations  are  repeated  (see  for  instance  Ref.  2) . 
One  starts  from 

-  (y+1)  d>  <b  +  d)  =0 
'  Yxyxx  Yyy 

(for  simplicity  we  have  set  a*  =  0) 
and  sets 

<PX  =  n  ,  <!>y  =  0  (33) 

<j>(x,y)  =  +  xn  +  y0  (34) 


Then 


<l>  = 


in 

'n  3x 


—  0  •— r 0  —  *t  n  +  X  1  a. 

"  9  3x  n  x  +  y 


in 

3x 


30 

3x 


,  _  .  3n  .  30  ^  ±  ,  30 

<P  =  -4>  +X7r-l+0  +  y-^— 

y  n  3y  0  3y  3y  3  3y 

Hence  with  Eq.  (33) 


0  -  H<-*„  +  *>  +  M'-*e  + 


0  =  +  x)  +  -(-4>  +  y) 

3y  n  3y  0 


and 


x  = 


<t> 

n 


y  =  % 


(35) 


23 


The  derivatives  <j>  ,  4>  and  4>  are  found  from  the 

xx  xy  yy 

identities 

n  =  $x(x(n,e)  ,y(r\,e) ) 

and 

e  =  (x ( n » 0 )  »y(n,e) ) 

One  obtains  by  differentiating  the  first  one  with  respect 
to  n  and  0 


1=4,  +  (j, 

yxx  3n  xy  3n 


0  ■  *xx  If  *  *xy  ft 


(36) 


(37) 


(38) 


(39) 


The  lowest  order  approximation  to  the  basic  flow  in  the  vicinity 
of  a  point  of  the  sonic  line  is  given  by 
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1 


<t>  (n,  6) 


R[ne  +  (l±i  n3  +  L_)  tgB) 


Hence 

x  =  <j>  =  R [ 0  +  n^tgB] 

n  *• 

y  =  $0  =  R[n  +  6tg6] 


(40) 


This  expression  satisfies 

$ 

A  _ _ Hi  =  -R _ _ 

XY  °  tg20R5(Y+l)n  -  R2 

Hence,  for  n  small 


in  accordance  with  Eq. 
There  one  has 


<J> 

xy 
(2)  . 


The  sonic  line  is  given  by  n 


0. 


x  =  R8 
y  =  R8tg6 

In  accordance  with  the  definition  following  Eq.  (2), 
one  obtains,  indeed, 


dy 


sonic 


/dx  =  tgB 


If  one  superimposes  to  4>  a  solution  $  then,  4>  and 

( 1 )  n 

4> '  are  the  x  and  y  components  of  the  displacement  vector. 

°  (1)  (1) 

We  consider  now  a  solution  $  for  which  evaluated 

at  the  sonic  line  (n  =  0)  has  a  jump  as  one  proceeds  from 

positive  to  negative  values  of  0  and  is  constant  otherwise. 

The  choice  of  the  signs  is  made  in  accordance  with  the 

behavior  of  the  flow  field  on  the  upper  side  of  a  profile. 

At  the  sonic  line  the  angle  0  decreases  as  one  proceeds  in 

the  downstream  direction  and  R  is  negative. 
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Particular  solutions  of  this  character  have  the  form 

(41) 


•<n,6>  =  n3/2f(S) 


where 


5  =  |(Y  +  D"1  “T 


(42) 


The  characteristics  of  Eq.  (39)  are  recognized  if  one  rewrites 
it  in  the  following  manner 


They  are  given  by 


=  ;<rH)1/V/2 


or 


0  +  (y+1)  (2/3)  r^2  =  const 


The  lines  £  =  1  are  the  characteristic  through  the  origin. 

One  solution  for  which  =  const  along  a  line  i;  =  const 
is  obviously  given  by 

<t>  =  (3/2)  (y+l)1/20 

which  leads  to 

f  =  C1/2  (43) 

but  it  is  not  suitable  for  the  present  purpose.  One  obtains, 
by  substituting  Eq.  (41)  into  Eq.  (39) 

CU-l)f"  +  [§  -  j)f’  +  Yjf  =  0  (44) 
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This  is  indeed  satisfied  by  Eq.  (43)  .  Eq.  (44)  is 
reduced  to  a  first  order  equation  by  means  of  Eq.  (43) . 


f  =  S1/2gU) 


5(5-Dg"  +  (y5  -  §)g'  =  o 


Hence 

2  _  1 

g"  _  2  3^  _  311  1 

g’  5(5~1)  "  ”  2  ?6  ?  -  1 


Then,  except  for  a  factor 


g '  =  r3/2(C  -  1)1/6 


g(5)  =  f  Z 
Z=i 


-3/2  ({.ui/Saj 


With  this  choice  of  the  constant  of  integration  one  obtains 
g  =  0,  g'  =  0  for  5  =  1.  This  gives  an  admissible  transition 
along  the  characteristic  5=1  (Fig.  8  in  the  second  quadrant) 
from  the  solution  <J>  ^  =  0 ,  to  the  expression  <i>  =  n3//2/5g(5) 
for  4>  as  well  as  its  first  derivative  vanish  along  this  line. 
The  sonic  line  is  given  by  n  =  0,  5  =  <».  One  has 


g  (°°) 


00 

/  z3/2  (Z-l)1/6dZ 
Z=l 


This  constant  can  be  evaluated  in  terms  of  tabulated  functions. 
One  obtains  by  setting 


Z  = 


5 


-1 


g  (°°)  = 


l 

/ 

o 


.1/3  -  1 


(1-5) 


7/6  -  1 


d5 


'  B<J-  5> 
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where  B  is  the  Eulerian  Beta-function.  Hence 


g(°°)  = 


r<7>  r<g> 

n|> 


One  has  for  £  small 


g(5_1)  =  /  r2/3(i-rj1/6d£ 


/  -  /  r2/3(i-U1/6<H 


r(i)  n|> 

r(|) 


-  3C 


1/3  _ 


r  (|) 


-  3r"1/3 


Hence  for  the  vicinity  of  the  line  n  =  0  ( £  =  °°) 


(1)  =  n3/2  .1/2  gU) 


3  -1/2  r4)T(J) 

|<Y+1)  1/2  -  3  y-g-  e 

F(I) 


34/3  2-l/3 (Y+1)-l/60l/3n 


Therefore,  for  n  =  0 


=  -34/3  2-1/3(y+1)-1/6  e1/3 


(1)  _  3 


=  t<y+D 


-1/2 


r<!> 


n|) 


The  displacement  in  the  y  direction  is  given  by  ^  .  In 
the  lowest  approximation  this  displacement  is  constant  at  the 
boundary  n  =  0  of  the  shaded  region  of  Fig.  8,  while  it  is 
zero  for  9  positive.  The  component  in  the  streamline  direction 
of  the  displacement  vector  is  then  given  by 
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(Y+l) 


(1)  _ 


=  -<{, 


(1) 


r{|) 


r<T> 


r(|> 


1/3  3 1/3  22/3  0l/3 


Now  we  had  in  the  basic  flow,  in  the  lowest  approximation, 
according  to  Eq.  (40) 


0 


x 

R 


Hence,  for  points  of  the  sonic  line  immediately  downstream  of 
the  point  where  the  jump  occurs 


dt  = 


-d 


n|) 


n  r (j)  r(|) 


( Y+l) 1/3  31/3  22/3  R"1/3  x1/3 


This  equation  displays  the  character  of  dt(x)  taken  along  the 

sonic  line  for  the  vicinity  of  a  point  where  d  (s)  has  a  jump. 

n  1/3 

Important  for  the  purpose  of  interpolation  is  the  factor  x  '  . 

Sometimes  it  may  not  be  necessary  to  determine  the  coefficient 
by  which  this  factor  is  multiplied.  The  qualitative  behavior 
of  some  quantities  which  are  needed  for  subsonic  computations 
is  shown  in  Fig.  9. 

The  singularity  generated  by  a  jump  of  dn(x)  will 
propagate  in  the  downstream  direction  along  the  characteristic 
which  starts  at  the  point  where  the  jump  occurs.  It  is 
reflected  at  the  surface  of  the  profile.  A  singularity  of 
the  same  character  will  occur  where  this  reflected  wave 
reaches  the  sonic  line.  This  wave  is  reflected  again.  These 
discussions  show  the  character  of  the  singularities  which 
the  kernel  K(s,t)  in  Eq.  (6)  will  have.  As  was  mentioned 
above,  the  numerical  evaluation  gives  K(s,t)  pointwise. 
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In  deriving  boundary  conditions  for  the  subsonic  region 
from  the  Eq.  (6) ,  one  must  discretize  the  relations  for 


Ad>  and  A <b  and  d  (s)  in  a  form  which  is  suitable  for  the 
Yx  y  n 


computation  of  the  subsonic  region.  This  means  that  one 
must  approximate  these  data  in  a  form  analogous  to  the 


expressions  for  A<t>^  and  Ad>y  which  appears  in  the  subsonic 


finite  element  formulation.  Details  will,  of  course,  depend 
upon  the  character  of  the  shape  functions  chosen  in  the 
subsonic  region. 

This  provides,  at  least  in  principle,  the  system  of 
formulae  needed  to  carry  out  the  procedure  sketched  in 
Section  II. 


1 


1 
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Figure  1.  Characteristics  Within  the  Field. 
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Figure  2.  Characteristic  Terminating  at  the  Sonic  Line. 


Finure  3.  Characteristic  Terminating  at  the  Profile  Contour. 
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Figure  4.  Characteristics  Connecting  Points  of  a  Basic  Field  (Points 
1,  2  and  3)  and  of  a  Perturbed  Field  (Points  1 '  ,  and  3'} 


Computation  of  the  Displacement  Vector  of  a  Point  of  the 
Sonic  line.  (Ax,,  Ay,  and  d  are  known,  d  is  to  be 
determined) 


Figure 
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Figure  6 


Figure  7 


Figure  8 


y.n 


Sketch  Illustrating  the  Relation  Between 


d6 

dn 


and 


d0 
dt ' 


de 

dy 


t 

r 


I  ^SONIC  LINE 

/ 

r7'r7PROFiLE  CONTOUR 


Starting  Point  (1)  for  the  Method  of  Characteristics . 


Hodograph  (ri,  8  plane).  The  Curves  c  =  1  are  the 
Characteristics  through  the  Origin.  The  Shaded  Part 
is  Affected  by  a  Jump  of  dn> 
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Figure  9 


Qualitative  Behavior  of  d  ,  A4>  and  A<J>  at  a  Point  of 
the  Sonic  Line  where  d  Jumps  y  from  ziro  to  a  Finite 
Value.  n 
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